This directory contains the python code that fits the rotational distortion constant to the spectroscopy from the main figure.

Its input is a hyperfine-containing model that has Dv set to 0, as the term proportional to Dv is diagonal in N we can add this
term on later. 

The output is to the python terminal so this should be run in such a way that this is visible.